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TRANSPORT  OF  WATER  IN  FROZEN  SOIL 
I.  Experimental  determination  of  soil-water 
diffusivity  under  isothermal  conditions 

Y.  Nakano,  A.R.  Tice,  J.L.  Oliphant  and  T.F.  Jenkins 


INTRODUCTION 

It  is  known  (Anderson  et  al.  1978)  that  water  in 
frozen  soil  generally  exists  in  three  phases:  vapor, 
liquid  water  and  ice.  The  liquid-phase  water  in  fro¬ 
zen  soil,  which  is  believed  to  be  mobile,  is  often  re¬ 
ferred  to  as  unfrozen  water  (Anderson  et  al.  1978). 
The  unfrozen  water  content  depends  on  the  soil  type 
and  temperature,  and  monotonicaDy  decreases  with 
decreasing  temperature.  It  is  known  that  some  clay 
contains  up  to  5%  unfrozen  water  by  weight  at 
-10°C.  Although  it  has  been  widely  assumed  that 
the  unfrozen  water  content  is  independent  of  the 
ice  content,  recent  work  (Tice  et  al.  1978)  using 
the  nuclear  magnetic  resonance  technique  indicates 
that  unfrozen  water  noticeably  increases  with  in¬ 
creasing  ice  content. 

It  has  been  recognized  that  accurate  knowledge 
of  water  movement  in  frozen  soils  is  very  important 
for  many  practical  applications.  For  example,  frost 
heave  is  known  to  be  caused  by  transport  of  water 
through  freezing  soils  (Loch  and  Miller  1975)  and 
the  damage  of  frost  heave  is  of  major  economic  con¬ 
cern  to  both  agriculturalists  and  engineers. 

Wiliams  and  Burt  (1974)  measured  the  hydraulic 
conductivity  of  frozen  soil  using  a  permeameter  in 
which  frozen  soil  was  positioned  between  two  com¬ 
partments  filled  with  an  aqueous  solution  of  lactose. 
However,  the  presence  of  significant  concentrations 
of  lactose  inevitably  introduces  a  deviation  in  the 
system  to  some  unknown  degree.  In  order  to  over¬ 
come  this  difficulty  Horiguchi  and  Miller  (1980) 
tatroduced  a  different  kind  of  permeameter,  in  which 
frozen  soil  was  positioned  between  two  compartments 


filled  with  pure  supercooled  water,  and  measured 
the  hydraulic  conductivity  of  frozen  soil  in  the  tem¬ 
perature  range  of  0°  to  -0.1 5°C.  The  lowest  temper¬ 
ature  attainable  in  this  permeameter  is  evidently  re¬ 
stricted  by  freezeup  of  the  supercooled  water. 

In  the  two  methods  of  permeability  measurement 
described  above,  the  lowest  attainable  temperature 
is  restricted  one  way  or  another.  In  this  paper,  we 
will  introduce  a  new  and  simple  method  that  could 
be  used  to  measure  the  rate  of  water  transport  in 
frozen  soil  at  any  temperature.  We  will  present  the 
theoretical  justification  of  this  new  method  and  some 
experimental  data  to  demonstrate  the  feasibility  of 
the  method. 

When  we  consider  water  transport  in  frozen  por¬ 
ous  media  under  isothermal  conditions,  it  is  reason¬ 
able  to  expect  that  this  problem  can  be  treated  in  a 
manner  similar  to  the  problem  of  unfrozen  soil  under 
the  condition  that  either  ice  is  absent  or  the  effects 
of  ice  on  water  transport  are  negligible.  Under  such 
conditions  water  transport  in  frozen  soil  corresponds 
to  that  in  unfrozen  soil  in  which  the  volumetric  water 
content  is  less  than  10%  or  so. 

Historically  in  soil  physics,  the  problem  of  water 
transport  in  soil  containing  a  small  amount  of  water 
has  been  subjected  to  intensive  investigations  (Swartz- 
endruber  1969).  If  we  assume  that  the  medium 
obeys  the  extended  Darcy's  law  for  unsaturated 
porous  flow  and  that  gravitational  effects  are  negli¬ 
gible,  the  flux  of  water  Fin  one  direction  under  an 
isothermal  condition  is  given  as 

F*-D(0)p&  (1) 


THEORY 


where  0  is  the  weight  of  water  per  weight  of  dry 
soil,  x  is  the  space  coordinate,  p  is  the  dry  density 
and  0(0)  is  the  soil-water  diffusivity  (we  will  use 
the  water  content  based  upon  weight  throughout 
this  work).  When  the  transport  of  water  in  the 
vapor  phase  is  included,  0(0)  is  given  (Philip  1957) 

3 

W-  E  A(«)  (2) 

i=l 

where  Dx ,  0a  and  03  are  the  diffusivities  of  ab¬ 
sorbed  liquid  water,  vapor,  and  pore  liquid  water, 
respectively.  Raats  (1975)  has  investigated  eq  2  by 
using  the  mixture  theory  of  continuum  mechanics 
and  has  shown  the  validity  of  eq  2  if  the  motion  of 
the  gaseous  phase  in  the  system  does  not  significantly 
affect  the  transport  of  water. 

In  frozen  soil  almost  all  pore  liquid  water  is  be¬ 
lieved  to  be  frozen  (Anderson  et  al.  1978)  and  we 
have 

D(9)  -0,(0)  +  0,(0).  (3) 

It  is  anticipated  that  the  transport  of  water  in  frozen 
soil  occurs  in  two  phases,  absorbed  liquid  water  and 
vapor,  and  that  an  overall  diffusivity  0(0)  is  given 
by  the  simple  sum  of  the  diffusivities  of  both  phases. 

When  eq  1  holds  true,  the  equation  of  continuity 
is  given  as 

If- £[«•>!]  «> 

where  t  is  time.  The  commonly  used  method  for 
the  experimental  measurement  of  0  utilizes  the 
mixed  initial  and  boundary  value  problem  of  eq  4 
with  initial  and  boundary  conditions  given  as 


0(x,O)»0, 

x*0 

«0j 

0<x 

(5«) 

0(O,/)-0t 

0<  r 

(5b) 

where  0,  and  02  are  constants  with  0,  >  03.  In  die 
reported  experiments  the  water  content  0,  is  kept 
constant  by  either  a  fritted  glass  bead  plate  (Nielson 
et  aL  1962)  or  a  diffusion  dumber  (Jackson  1965). 

in  this  work  we  will  introduce  the  method  of  ex¬ 
perimental  measurement  by  using  the  initial  value 
problem  of  eq  4  in  order  to  avoid  the  inconvenience 
of  maintaining  the  water  content  constant  at  the 
boundary. 


Our  experiment  consists  of  two  long  columns  of 
seal  with  the  same  size  and  dry  density.  One  of  them 
is  uniformly  dry  with  the  negligibly  small  water  con¬ 
tent  0  =  0  while  the  other  is  uniformly  wet  with 
0  =  0O  >  0.  The  wet  column  is  allowed  to  contain 
ice  only  for  the  case,  if  it  exists,  in  which  the  effects 
of  ice  on  water  transport  are  negligible.  At  time 
t  =  0  we  connect  the  two  columns  to  make  a  single 
column  from  which  no  water  escapes.  While  we 
maintain  the  column  at  the  specified  constant  tem¬ 
perature,  water  is  transported  from  the  wet  part  to 
the  dry  part.  This  experiment  can  be  described  in 
mathematical  terms  by  the  following  initial  value 
problem  of  eq  4  on  the  set,  S  =  {(-“,  °°)x(0,  T); 
T>0\: 

HrH  <“> 

0(x,O)  =  0o  x  <  0 

=  0  x  >  0 .  (6b) 

Oleinik  et  al.  (1958)  investigated  the  general  in¬ 
itial  value  problem  of  eq  6a  and  have  shown  the  ex¬ 
istence  and  uniqueness  of  the  solution  under  the 
conditions  that  0(x,  0)  satisfies  the  Lipschitz  con¬ 
dition  for  -*>  <  x  <  *»,  that  0(0)  has  continuous 
derivatives  of  the  fourth  order,  and  that  0(0)  is 
bounded  for  bounded  0.  It  should  be  mentioned 
that  the  solution  defined  (Oleinik  et  al.  1958)  may 
not  be  continuously  differentiable  and  is  often  re¬ 
ferred  to  as  a  weak  (generalized)  solution.  Unfortun¬ 
ately,  the  existence  and  uniqueness  theorem  (Oleinik 
et  al.  1958)  is  not  directly  applicable  to  the  problem 
of  eq  6a  and  b  because  the  initial  condition  of  eq  6b 
obviously  does  not  satisfy  the  Lipschitz  condition. 

Recently,  Nakano  (in  press)  investigated  a  special 
case  of  the  problem  of  eq  6a  and  b  in  which  0(0) 
is  given  as 

0(0)-0o-0o0*  (7) 

where  0O  and  0  are  positive  numbers.  Equation  6a 
can  be  transformed  into  an  ordinary  differential  equa¬ 
tion  by  introducing  the  similarity  variable  q  *  x  [0O 
(/  +  f0)l  where  f0  >  0  is  a  constant.  The  similarity 
solution  0(x,  t)  mAv)  should  satisfy  the  equation 

\f*f]  ’  +  (l/2)t/'-0  -»><  n<~  (8a) 
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where  primes  denote  differentiation  with  respect  to 
r>.  The  solution  /(17)  is  required  to  satisfy  the  condi¬ 
tions 

/(-«)- «o  /(°°)  =  0 .  (8b) 

Nakano  (in  press)  has  shown  that  there  exists  at 
least  one  continuous  solution  of  eq  8a  and  b.  This 
solution  is  strictly  monotone-decreasing  and  has  the 
property  given  as 

/(ij)  >0  for  a  >  T) 

f(rj)  =  0  for  a  <  rj  (9) 

where  a  is  a  positive  number.  We  will  call  the  inter¬ 
face  between  the  part  where  6  >  0  and  the  part 
where  $  ■  0  a  wetting  front.  It  follows  from  eq  9 
that  the  location  of  the  wetting  front  X(t)  is  given 

as 

*W  =  0/>*(f+fo)*.  (10) 

For  convenience’s  sake  we  define  the  functions, 
/,(d)  and/,  fri),  as 

/,(»»)  =Av)  if  t0  =0 

/2(d)=/(d)  if  to>0-  (11) 

It  has  been  shown  (Nakano,  in  press)  that  the  solution 
/,(»f)  is  not  unique  and  contains  one  unspecified  pa¬ 
rameter.  Choosing  a  to  be  this  unspecified  parameter, 
we  will  write /,(»?)  as/, (ij;o).  This  one-parameter 
family  of  solutions/,  (ij;  a)  satisfies  the  initial  con¬ 
dition  equation  (eq  6b)  and  also  provides  solutions 
to  the  problem  of  eq  6a  and  b,  implying  that  the 
problem  of  eq  6a  and  b  has  many  solutions.  The 
solution /3(t;)  does  not  satisfy  the  initial  condition 
equation  (eq  6b).  Instead,  the  initial  data  of  f2  (i?) 
have  the  property  given  as 

/2(d)>  0  17  <  at* 

/2(d)«0  V>ot#.  (12) 

Unlike  /,  (17;  o),/3(p)  is  shown  to  be  unique. 

Studying  the  asymptotic  behavior  of  die  solution 
to  the  wide  class  of  initial  value  problems  of  eq  6a, 
Nakano  (in  press)  has  shown  that  the  solutions  of 
eq  6a  and  b,/,(d;  o),  asymptotically  converge  upon 
/3(d)  as  time  tends  to  infinity  with  the  proper  choice 
of  a.  The  proper  value  of  0  cannot  be  determined 


from  the  initial  condition  equation  (eq  6b),  but  from 
the  experimental  data  as  described  below. 

The  properties  of  the  solution  of  eq  8a  and  b, 

/( 17),  were  investigated  (Nakano,  in  press)  and  it  has 
been  shown  that  /2(d)  is  not  continuously  differen¬ 
tiable  because  the  transition  between  the  part/>  0 
and  the  part /=  0  is  not  smooth.  In  left  neighbor¬ 
hoods  of  the  wetting  front  (o  -  e  <  1}  <  o;  6  >  0)/(tj) 
is  found  to  behave  as 

m=(omv/e(o-riytfi .  03) 

The  derivatives  of /(1?)  are  given  as 

/'(d)  =  -r 1  (o0/2)I,e  (o  -  V)°  -  M  (14a) 

An) =0  -d)r2(o0/2)*/^(o-T,)o-^ 

(14b) 

Ad)  =  -0  -0X1  -2/3)r3(ad/2)«/U(a-T7)0 

(14c) 

etc.  From  the  above  equations  we  find  the  following. 
If  P  >  1 ,  Ad)  tends  to  minus  infinity  as  tj  tends  to  a 
and  Ad)  <  0.  This  implies  that  the  wetting  front  is 
very  steep  and  convex  upwards.  This  type  of  wetting 
front  was  reported  (Nielson  et  al.  1962).  If  1  >  0  >  0, 
Ad)  remains  finite  and  Ad)  >  0.  This  type  of 
wetting  front  was  observed  (Jackson  1 964).  It  is 
clear  that  the  transition  between  the  part  /  >  0  and 
the  part  /*  0  becomes  smoother  as  0  decreases.  When 
0  =  0,  it  is  known  that  there  is  a  continuously  differ¬ 
entiable  solution  given  as 

0(d)  *  (0o/2)  (1  “  erf  (d/2)]  (15) 

where  erf  is  the  error  function  defined  as 

1  * 

erf(x)=  — —  f  exp  (-sJ)  ds.  (16) 

Vff  0 

It  follows  from  eq  1 5  that  9  tends  to  zero  as  17 
approaches  plus  infinity.  Hence,  according  to  our 
definition,  there  appears  no  wetting  front  in  this 
solution.  Equation  15  may  be  written  as 

0O  “  0(d)  =  (0o/2)  erfc  (-d/2)  (17) 

where  erfc  is  the  complemental  error  function  de¬ 
fined  as 
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erfc  (x)  *  1  -erf(x) .  (18) 

When  t;  approaches  minus  infinity,  the  solution  /(r?) 
has  beat  shown  (Nakano,  in  press)  to  behave  as 

«o  "/(n)  »  («0/2)  «fc  l-u/(2flf)l  .  (19) 

The  behavior  of f(rj)  resembles  d(tj)  near  tj  =  -os. 

When  D(B)  is  given  as  in  eq  7,  we  would  measure 
the  profile  of  water  content  with  time  in  our  experi¬ 
ment.  It  is  expected  that  after  some  time  t-T,  the 
profile  will  become  stable  and  vary  little  with  time 
afterward.  We  take  the  profile  at  t  *  T  as  an  experi¬ 
mental  asymptotic  curve  f*(n).  By  comparing/^!?) 
with  the  theoretical  profile  /t(i);  a)  we  can,  in  prin¬ 
ciple,  determine  the  values  of  a,  D0  and  0.  For  in¬ 
stance,  from  the  measured  locations  of  a  wetting 
front,  X(T),  we  find 

aD$^X(T)r-».  (20) 

In  the  next  step  we  seek  approximate  values  of 
D0  and  0  using  eq  13.  Equation  13  may  be  written 

as 


Z*ln/*a  +  0->  Y 

(21) 

««0-‘b(o0D5»/2) 

(22a) 

y-aDj-xr-*. 

(22b) 

Using  several  experimental  points  (x,f)  near  the 
wetting  front,  we  plot  Z  vs  Y.  The  graphs,  Z)is 
expected  to  be  approximately  linear.  Thus  we  can 
determine  the  approximate  values  of  D0  and  0. 

In  the  final  step  we  seek  the  values  of  D0  and  0 
that  minimize  the  difference  between  /•(»?)  and 
/,(tj;  o)  by  trial  and  error.  The  theoretical  profile 
/j(t};  a)  can  be  obtained  by  the  numerical  integration 
of  eq  8a.  We  will  describe  below  one  of  several  pos¬ 
sible  methods  of  integration  that  we  use  in  this  work. 

For  convenience’s  sake  we  introduce  two  new 
variables  defined  as 

s  ■  n/o  (23a) 

yffoW.  (23b) 

Using  eq  23a  and  b,  we  reduce  eq  8a,  13  and  14a  to 
0*/y  +  Vxy'  ■  0  -~  <  i  <  1  (24) 


y  =  (0/2)«/0(i  -*)«/<» 

(25) 

y  =  -0-'  (0/2)1#  (1  -j)<‘-0>/0 

(26) 

where  primes  denote  differentiation  with  respect 
to  s.  In  order  to  reduce  eq  24  to  a  system  of  first- 
order  equations,  we  introduce  another  dependent 
variable  defined  as 

z  =  -yty  . 

(26) 

Now,  eq  24  is  reduced  to 

y'  =  ~y-Bz  =  F(s,y,z) 

(27a) 

z'  =  -  Visy-tz  =  G  (s,y,  z) 

(27b) 

with  conditions 

y(l)  =  z(l)  =  0. 

(28) 

We  will  simultaneously  integrate  eq  27a  and  b  by 
the  Runge-Kutta  method  from  s  =  1  to  any  desired 
value  of  s  <  1  using  eq  25  and  26  for  evaluating  F 
and  G  at  the  first  step  of  integration.  Since  y  varies 
most  rapidly  near  the  wetting  front  s  -  1  and  the  rate 
of  change  decreases  with  decreasing  s,  it  is  desirable 
to  use  variable  increments.  Particularly,  when  0  is 
greater  than  1 ,  the  selection  of  the  first  step  incre¬ 
ment  should  be  carefully  examined  in  order  to  attain 
the  required  accuracy  of  computation.  In  the  present 
work  the  value  of  0  is  much  less  than  1  and  we  com¬ 
puted  the  profile  of  water  content  with  a  relative 
error  of  less  than  1%  using  the  initial  increment  of 
10-*. 


EXPERIMENT 

The  soil  selected  for  this  study,  called  Morin  clay, 
is  a  marine-deposited  clay  obtained  from  the  Morin 
brickyard,  Auburn,  Maine.  This  area  was  uplifted 
and  the  original  solutes  leached  by  natural  seepage. 
The  specific  surface  area  is  60  m3/g  soil,  which  may 
be  considered  typical  for  clay  of  this  texture.  We 
present  the  measured  unfrozen  water  content  in  Morin 
clay  with  21%  total  water  content  at  various  temper¬ 
atures  together  with  the  fitted  curve  in  Figure  1 . 
Figure  1  shows  that  about  10%  of  the  water  remains 
unfrozen  at  -1 .0°C,  the  temperature  selected  for  this 
study. 

We  used  10-mL  graduated  plastic  pipettes  to  en¬ 
close  the  soil  columns.  The  dispensing  end  of  each 
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Figure  1.  Unfrozen  water  content  of  Morin  day  with  21% 
total  water  content  versus  temperature. 


pipette  was  cut  off  and  sealed  by  a  stopper.  The 
dimensions  of  each  soil  column  were  0.80  cm  in 
diameter  by  20.33  cm  in  length.  Each  experiment 
of  this  study  required  the  preparation  of  two  soil 
columns  with  the  same  specified  soil  dry  density. 

One  of  them  was  uniformly  dry  with  a  negligibly 
small  water  content  while  the  other  was  uniformly 
wet. 

We  chose  the  average  soil  dry  density  to  be  1.35 
g/cm3  and  the  water  content  of  the  wet  column  to 
be  8.00%  in  this  investigation.  Since  all  experiments 
were  conducted  under  the  isothermal  condition  of 
-1.0°C,  it  is  dear  from  Figure  1  that  no  ice  was 
present  in  this  study.  For  the  preparation  of  wet 
columns,  we  added  16  g  of  distilled  wateT  to  200  g 
of  oven-dry  Morin  day.  Mixing  thoroughly  the 
water  and  soil,  we  allowed  the  mixture  to  set  for  at 
least  3  days  to  attain  moisture  equilibration. 

The  uniform  pecking  of  soil  into  a  pipette  is  one 
of  the  most  critical  operations  in  this  experimental 
work.  Despite  frequent  weighing  of  die  soil-filled 
pipettes,  we  had  to  discard  some  of  them  due  to 
lack  of  uniformity  in  terms  of  dry  density.  For  the 
preparation  of  dry  columns,  oven-dry  Morin  cky  was 
packed  into  a  pipette  with  slightly  greater  compac- 
tive  force.  During  the  packing,  Morin  clay  gained 
a  small  amount  of  water  from  the  air,  and  the  water 
content  of  dry  columns  avenged  about  0.4%  at  the 


end  of  packing.  Immediately  after  the  packing  was 
finished,  we  sealed  the  open  end  with  a  sleave  type 
of  rubber  stopper  and  allowed  the  sample  to  set  for 
at  least  3  days  for  moisture  reequilibration.  Wet 
columns  were  frozen  down  to  -16.0°C  and  gradually 
warmed  up  to  -1 ,00°C. 

An  experiment  began  by  joining  two  pipettes,  one 
with  dry  soil  and  one  with  wet  soil,  which  were  both 
set  at  -1.0°C.  Upon  removing  the  sleave-type  stop¬ 
pers  from  these  two  pipettes,  we  joined  them  with 
surgical  rubber  tubing.  Both  surfaces  of  contact 
were  made  smooth  in  order  to  avoid  contact  resis¬ 
tance.  Immediately  after  the  two  pipettes  were 
joined,  they  were  placed  in  a  constant  temperature 


Table  1.  Mean  (g/cm1)  and  standard  deviation  of 
dry  density. 


Exp.  series 

Time  duration 
(dart 

Mean 

Standard 

deviation 

1 

2 

1.35 

0.102 

4 

1.36 

0.120 

8 

1.3$ 

0.111 

16 

1.35 

0.131 

I 

4 

1.36 

0.122 

8 

1.34 

0.123 

16 

1.35 

0.116 

32 

1.35 

0.115 
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bath  containing  an  ethylene  glycol-water  mixture. 
The  bath  temperature  was  set  at  -1 .0°C  and  was 
maintained  to  within  ±  0.03°C  using  a  Bayley  pro¬ 
portional  temperature  controller. 

After  a  specified  time  passed,  the  jointed  pipette 
was  removed  from  the  constant  temperature  bath 
and  quickly  sectioned  into  a  total  of  64  equal  thin 
column  segments.  Each  segment  was  placed  in  a 
glass  weighing  bottle  and  oven-dried  overnight.  The 
water  content  and  die  dry  density  of  each  segment 
were  determined  gravimetrically. 

The  time  durations  were  2, 4,  8  and  16  days  for 
the  first  series  of  experiments.  In  order  to  examine 
the  reproducibility  of  the  experiments,  we  conducted 
the  second  series  of  experiments  where  the  experi¬ 
mental  times  were  set  4, 8, 16  and  32  days. 

The  uniformity  of  the  soil  columns  was  examined 
by  using  the  measured  dry  density  of  segments  after 
each  experiment.  The  mean  and  the  standard  devia¬ 
tion  of  dry  density  calculated  are  given  in  Table  1 . 


RESULTS  AND  DISCUSSION 

We  present  the  experimental  results  in  Figure  2 
where  the  water  content  0  is  plotted  versus  the  dis¬ 
tance  x.  As  described  above  we  have  two  sets  of 
data  at  f  =  4,8  and  16  days.  These  two  sets  of  pro¬ 
files  in  Figure  2  indicate  that  it  is  difficult  to  accur¬ 
ately  reproduce  the  same  profile,  and  that  the  dif¬ 
ference  *•  'tween  the  two  profiles  amounts  to  about 
0.5  in  terms  of  percentage  of  water  content,  which 
corresponds  to  about  10%  relative  error  in  the  wetter 
end  of  a  soil  column. 

According  to  the  results  of  the  mathematical 
analysis  described  in  the  previous  section,  experi¬ 
mental  profiles  are  expected  to  converge  upon  an 
asymptotic  profile  as  time  increases  if  D(0)  is  given 
by  eq  7.  In  order  to  examine  the  asymptotic  be¬ 
havior  of  experimental  profiles,  in  Figure  3  we  plot¬ 
ted  the  water  content  6  versus  £  given  as 

(29) 

In  Figure  3  we  used  the  average  values  of  two  sets 
of  data  at  t  *  4, 8  and  16  days.  This  figure  indicates 
that  the  experimental  profiles  coincide  with  each 
other  reasonably  well  after  8  days.  We  conclude  that 
experimental  profiles  tend  to  stabilize  as  if  they  were 
converging  upon  an  asymptotic  profile  as  time  in¬ 
creases. 

In  the  next  step  we  will  examine  how  well  eq  7 
expresses  the  actual  behavior  of  D(0)  of  Morin  clay. 
Since  die  initial  water  content  of  dry  columns  of 


x(cm) 


Figure  2.  Water  content  profiles:  (1)1  =  2 
days,  (2)1  =  4  days,  (3)  t  =  8  days,  (4)t  =  16 
days,  (5)  t  =  32  days. 

soil  is  not  negligibly  small  in  our  experiments,  we 
introduce  a  new  variable  u  given  as 

u  =  0-0c  (30) 

where  0C  is  the  initial  water  content  of  dry  columns. 
Using  eq  30  we  reduce  eq  7  to 

D(d)=D0(e-ecy.  (3i) 

It  is  easy  to  find  that  the  results  of  the  mathematical 
analysis  described  in  the  previous  section  hold  true 
in  terms  of  the  new  variable  u  if  D(0)  is  given  by 
eq  31. 

Using  the  least-squares  method  we  determined 
the  parameters  in  eq  31  for  the  averaged  experimental 
profiles  at  t  =  8  days  as 

D(6)  =  1.21  (D  -  3.8  x  10-5)0  271  (32) 

where  D(0)  is  in  cm2 /day  and  0  is  in  g  water/g  dry 
soil.  The  curve  in  Figure  3  is  the  computed  profile 


Figure  3.  Water  content  versus  £. 


Figure  4.  Soil-water  diffusivity  versus  water  content  (Adel- 
anto  and  Bachappa  loams  from  Jackson  ( 1965]  ). 


using  eq  32  and  the  Runge-Kutta  method  described 
earlier.  The  computed  profile  appears  to  be  a  rea¬ 
sonably  good  approximation,  but  is  not  quite  a  close 
fit  for  the  experimental  profiles. 

In  order  to  examine  the  reason  why  eq  31  does 
not  provide  a  close  fit,  we  will  evaluate  D(fi)  directly 
from  the  experimental  profile  at  t  =  8  days  under 
an  assumption  described  below.  Using  the  similarity 
variable  £,  we  reduce  eq  da  and  b  to 

[D(f)f T+%|/*-0  -«<£<«  (33a) 

/(-«)  =  «o  /(°°)  =  0 .  (33b) 

We  assume  that  the  solution  of  eq  6a  and  b  asymp¬ 
totically  converges  upon  that  of  eq  33a  and  b  re¬ 
gardless  of  the  functional  form  of  D(0 )  if  D(0)  >  0 
is  a  smooth  function  of  0.  When  this  assumption 
holds  true,  by  integratkig  eq  33a  we  obtain 


om i  -[mo-p+  j  H/[2_i] 

(34) 

where  £0  is  the  value  of  {  at  the  wetting  front.  Using 
the  experimental  profile  at  t  =  8  days,  we  have  ap¬ 
proximately  evaluated  the  right-hand  side  of  eq  34 
to  obtain  D(0).  We  have  plotted  the  calculated 
values  of  D(0)  versus  6  in  Figure  4  where  the  curve 
is  eq  32.  In  spite  of  inherent  inaccuracies  in  evalu¬ 
ating  derivatives  from  experimental  data,  it  appears 
that  there  is  a  slight  peak  in  D(0)  near  0  =  1%.  Due 
to  uncertainty  as  to  the  validity  of  this  result,  we 
can  only  speculate  that  eq  31  does  not  provide  a 
close  fit  for  experimental  profiles  because  the  actual 
D(0)  may  not  be  a  monotonically  increasing  function 
offl. 
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It  seems  that  no  one  has  reported  the  soil-water 
diffusivity  of  frozen  soils.  Jackson  (1965)  measured 
the  diffusivity  of  Adelanto  loam  and  Pachappa  loam 
at  temperatures  of  6.5°C  and  6.0°C,  respectively. 

We  plotted  Jackson’s  data  in  Figure  4.  It  is  known 
that  D(0)  decreases  as  temperature  decreases. 
Assuming  that  Jackson’s  theoretical  prediction 
based  upon  the  concept  of  activation  energy  can  be 
extended  to  -10°C,  we  roughly  estimate  that  the 
diffusivity  at  -1 .0°C  would  be  three-quarters  to 
one-half  of  that  at  6.0°C.  This  brings  Jackson’s 
and  the  present  data  closer  together,  and  our  meas¬ 
ured  values  of  the  soil-water  diffusivity  are  compar¬ 
able  to  Jackson’s  data. 


CONCLUDING  REMARKS 

We  have  introduced  a  new  experimental  method 
for  measuring  the  soil-water  diffusivity  of  frozen 
soil  under  isothermal  conditions.  We  have  presented 
the  theoretical  justification  of  the  method  and  dem¬ 
onstrated  the  feasibility  of  the  method  with  an  ex¬ 
periment  using  Morin  clay. 

The  theoretical  justification  is  made  under  the 
assumption  that  the  soil-water  diffusivity  is  propor¬ 
tional  to  some  power  of  water  content.  We  have 
found  that  this  functional  form  of  the  soil-water 
diffusivity  does  not  provide  a  close  fit  for  the  exper¬ 
imental  data  obtained.  Thus  we  should  consider 
the  measured  diffusivity  approximate.  The  approx¬ 
imation  was  needed  in  analyzing  experimental  data 
because  the  mathematical  theorems  currently  avail¬ 
able  are  limited.  Mathematical  theorems  for  the 
more  general  case  of  a  diffusivity  function  are  re¬ 
quired  to  improve  the  accuracy  in  data  analysis. 


LITERATURE  CITED 

Anderson,  D.M.,  R.  Pusch  and  E.  Fenner  (1978) 
Physical  and  thermal  properties  of  frozen  ground. 
Geotechnical  Engineering  for  Cold  Regions  (O.B. 
Andersland  and  D.M.  Anderson,  Eds.).  New  York: 
McGraw-Hill,  pp  37. 

Hotiguchi,  K.  and  R.D.  Miller  (1980)  Experimental 


studies  with  frozen  soil  in  an  ‘ice  sandwich’  perme- 
ameter.  Cold  Regions  Science  and  Technology, 

3:  177. 

Jackson,  R.D.  (1964)  Water  vapor  diffusion  in  rela¬ 
tively  dry  soil:  I.  Theoretical  considerations  and 
sorption  experiments.  Soil  Sciences  Society  of 
American  Proceedings,  28:  172. 

Jackson,  R.D.  (1965)  Water  vapor  diffusion  in  rela¬ 
tively  dry  soil:  IV.  Temperature  and  pressure  ef¬ 
fects  on  sorption  diffusion  coefficients.  Soils  Science 
Society  of  American  Proceedings,  29:  144. 

Loch,  JJ*.G.  and  R.D.  Miller  (1975)  Tests  of  the 
concept  of  secondary  frost  heaving.  Soils  Science 
Society  of  America"  Proceedings,  39:  1036. 

Nakano,  Y.  (in  press)  Asymptotic  behavior  of  solu¬ 
tions  of  the  problem  of  wetting  fronts  in  one-dimen¬ 
sional,  horizontal  and  infinite  porous  media.  Ad¬ 
vances  in  Water  Resources. 

Nielson,  D.R.,  J.W.  Biggar  and  J.M.  Davidson  (1962) 
Experimental  considerations  of  diffusion  analysis 
in  unsaturated  flow  problems.  Soils  Science  Society 
of  American  Proceedings,  26:  107. 

Oleinflc,  O.A.,  A.S.  Kalaahnlcov  and  C.  Yui-Lih  (1958) 
The  Cauchy  problem  and  boundary  problems  for 
equations  of  the  type  of  non-stationary  infiltration. 
Izv.  Akad.  Nauk.  USSR  Ser.  Mat,  22:  667. 

Philip,  J.R.  (1 957)  The  theory  of  infiltration :  1 .  The 
infiltration  equation  and  its  solution.  Soil  Sciences, 
83:  345. 

Raata,  P.A.C.  (1975)  Transformations  of  fluxes  and 
forces  describing  the  simultaneous  transport  of  water 
and  heat  in  unsaturated  porous  media.  Water  Re¬ 
sources  Research,  1 1 :  938 
Swartzendruber,  D.  (1969)  The  flow  of  water  in 
unsaturated  soils.  In  Flow  Through  Porous  Media 
(J.M.  DeWiest,  Ed.).  New  York:  Academic  Press, 
pp.  215. 

Tice,  A.R.,  C.M.  Burrows  and  D.M.  Anderson  (1978) 
Phase  composition  measurements  on  soils  at  very 
high  water  contents  by  the  pulsed  nuclear  magnetic 
resonance  technique.  Transportation  Research  Rec¬ 
ord,  No.  675,  pp.  11. 

Williams,  P.J.  and  T.P.  Burt  (1964)  Measurement  of 
hydraulic  conductivity  of  frozen  soils.  Canadian 
Geotechnical  Journal,  1 1 :  647. 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Nakano,  Y. 

Transport  of  water  in  frozen  soil.  I.  Experimental 
determination  of  soil-water  diffusivity  under  iso¬ 
thermal  conditions  /  by  Y.  Nakano,  A.R.  Tice,  J.L. 
Oliphant  and  T.F.  Jenkins.  Hanover,  N.H.:  U.S.  Cold 
Regions  Research  and  Engineering  Laboratory;  Spring- 
field,  Va.:  available  from  National  Technical  Infor¬ 
mation  Service,  1983. 

ii,  13  p.,  illus.;  28  cm.  (  CRREL  Report  83-22.  ) 

Bibliography:  p.  8. 

1.  Experimental  data.  2.  Frost  heave.  3.  Frozen 
soil.  4.  Soil  dynamics.  5.  Soil  water.  I.  Tice,  A.R. 
II.  Oliphant,  J.L,  III.  Jenkins,  T.F.  IV.  United 
States.  Army.  Corps  of  Engineers.  V.  Cold  Regions  Re¬ 
search  and  Engineering  Laboratory,  Hanover,  N.H. 

VI.  Series:  CRREL  Report  83-22. 


END 

FILMED 

1-84 


DTIC 


